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The transport properties of a Bose-Einstein condensate in a ID incommensurate bichromatic lat- 
tice are investigated both theoretically and experimentally. We observe a blockage of the center of 
mass motion with low atom number, and a return of motion when the atom number is increased. 
Solutions of the Gross-Pitaevskii equation show how the localization due to the quasi-disorder in- 
troduced by the incommensurate bichromatic lattice is affected by the interactions. 
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The intrinsic perfection of lattices made from a stand- 
ing wave of light surprisingly makes them an excellent 
candidate for the investigation of disorder in atomic sys- 
tems. Free from uncontrollable or undesirable defects, 
precise disorder can be added sirnply in the form of addi- 
tional optical lattices 
speckle potential @, [v 
cal lattices with a Bose-Einstein condensate (BEG) offers 
the stimulating complexity of interactions in a setting 
unimpeded by large thermal fluctuations. The signifi- 
cant observation of the Mott-Insulator phase was realized 
utilizing a strongly interacting BEG produced in a three- 
dimensional ordered crystal of light [Io| . Extending this 
work to the strongly disordered regime with the inclusion 
of quasi-disorder from a bichromatic lattice, recently led 
to initial experimental evidence of a Bose-glass phase [ll| . 
Another open question remains as to the intermediate be- 
havior between the non-interacting disordered Anderson 
localized phase and the strongly interacting Bose-glass 
phase [l^ . 

To begin to address this cross-over regime we have car- 
ried out investigations of the transport properties of a 
condensate in a quasi-periodic lattice. The very com- 
plexity provided by the interplay of disorder and inter- 
actions that makes disordered condensates a stimulating 
topic, also can introduce instability inherent in non-linear 
transport In fact we find that dynamical instability 
does occur for transport in a quasi-periodic lattice when 
the center of mass motion reaches a critical velocity, and 
has a non-trivial dependence on the interaction strength. 

Anderson's seminal paper in 1958 showed that the 
wave function of a particle placed in a lattice with disor- 
dered on-site energies remains localized when the range 
of the on-site energies is sufficiently large compared to 
the hopping energy between neighboring sites [ij. An- 
derson's hopping model can be approximated using a 
quasi-periodic bichromatic lattice, obtained by superim- 
posing a primary optical lattice with a weak secondary 
lattice of a different incommensurate wavelength. The 
secondary lattice serves to break the discrete transla- 
tional invariance of the system, thus favoring localization 
of the wave functions, however the effect of the quasi or- 



der may be important depending on the exact parameters 
of the bichromatic lattice [3, 0] ■ 

In our system the ID incommensurate bichromatic lat- 
tice is produced combining the primary optical lattice 
derived from a Titanium:Sapphire laser operating at a 
wavelength Ai = 830.7(1) nm with a secondary lattice 
obtained from a fiber-amplified diode laser emitting at 
A2 = 1076.8(1) nm. Our ^^Rb BEG is produced in a 
loffe-Pritchard magnetic trap, elongated in the direction 
of the bichromatic lattice. The trapping frequencies are 
cja; = 27r X 8.7 Hz axially and oj± = 27r x 90 Hz radially. 
The BEG can be produced in the range of ~ 1.5 x 10^ - 
2 x 10^ atoms. The resulting potential along the lattice 
axis can be expressed as 

V{x) = siEmsm'^{kix)-]-S2EB,2sin'^{k2x)-\-^ujlx'^ (1) 

where si and S2 measure the height of the lattice po- 
tentials in units of the respective recoil energies Em — 
/iV(2mAf) ~ /i X 3.33 kHz and Er2 = /iV(2mAi) ~ 
h X 1.98 kHz, ki and k2 are the wave numbers of the two 
lasers, h is the Planck constant and m the mass of ^^Rb. 

The possibility of localization with our bichromatic lat- 
tice in the absence of interactions is demonstrated by nu- 
merical diagonalization of the stationary ID Schrodinger 
equation using the potential defined by Eq. (1). A strong 
primary lattice is chosen, si = 10, siEm/h — 33 kHz, 
and a perturbing secondary lattice of maximum height 
S2 = 2, S2Eji2/h = 4 kHz. The ground state resulting 
from the bichromatic lattice is shown in Fig. 1(a) and 
is contrasted with the ground state of a pure random 
case in (b). The random potential is simulated using 
only the primary lattice, si = 10, with additional ran- 
dom on-site energies in a box distribution in the range [0, 
A], where A/h < 4 kHz. The amount of the disorder is 
given by either the height of the secondary lattice in (a) 
or by the maximum on-site energy in the random case 
(b), and is denoted by A in both (a) and (b). We de- 
fine J as the tunnelling in the primary lattice. The thin 
line is the ground state with only the primary lattice, 
showing the total length of the system in the harmonic 
trap. The localized states may have a different center 
from the harmonic trap depending on the disorder real- 
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FIG. 1: Non-interacting density profiles in log scale of the 
ground state with increasing disorder in (a) a bichromatic 
lattice and (b) a lattice with random on-site energies. The 
thin line in each graph represents the ground state with only 
the primary lattice. The last graph in each column shows the 
on-site energies in the respective potentials with the amount 
of disorder shown by A/J ~ 26. In all cases the height of 
the primary lattice is si = 10 giving a tunnelling energy of 
J/h = 75 Hz. 



ization. The last graph in each column is an enlarge- 
ment of the actual potentials. The quasi-periodic system 
mimics true disorder to a certain extent, showing local- 
ized states characterized by an exponential decay in the 
envelope of the density moving away from the localiza- 
tion center, oc exp(— |a; — xq\/1) where I is the 
localization length. The exponential localization occurs 
only above a threshold level at approximately A/J ~ 6 
in the bichromatic case. This is in clear contrast with 
random disorder where in a one-dimensional infinite sys- 
tem the localization behavior persists for any infinitesi- 
mal amount of disorder. 

In the presence of weak interactions one expects lo- 
calization effects to persist. In fact an Anderson glass 
phase has been identified when increasing interactions 
from an Anderson localized phase to the Bose-glass state, 
using the superfluid density as the order parameter for 
the phase transition. Unlike the Bose-glass phase, where 
both interactions and disorder cooperate to block diffu- 
sion, in the Anderson glass phase the interactions tend 
to delocalize the atoms and the disorder must compete 
against the interactions [13, [3] ■ 

Such behavior can be seen in Fig. 2, showing the 
ground states in our incommensurate bichromatic lat- 
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FIG. 2: Density profile in log scale of the ground state with 
interactions for (a)-(c) increasing atom number and fixed 
A/J = 6.5, and for (d)-(f) increasing A/J and fixed atom 
number iV = 1.5 x 10*. In aU cases the height of the primary 
lattice is si = 10 giving a tunnelling energy of J/h = 75 Hz. 



tice with interactions calculated by means of a ID ef- 
fective Gross-Pitaevskii equation, namely the nonpoly- 
nomial Schrodinger equation (NPSE) [iGj . 

In the presence of interactions the strongly localized 
ground state of (a) transforms into a state with multiple 
peaks with partially overlapping tails (b). Upon increas- 
ing the interactions, the overlap between these peaks in- 
creases until the state eventually becomes extended (c). 
These results show similar behavior toprevious simula- 
tions of a BEG in a three-color lattice . Even with the 
increased interactions, a state with multiple peaks can 
be recovered by increasing disorder (d)-(f). The cross- 
over behavior between an extended superfluid state and 
a possible Anderson glass state is very difficult to quan- 
tify in this picture. We take a pragmatic approach and 
characterize the state of the system by investigating the 
transport properties. 

In the experiment, to set in motion our harmonically- 
trapped condensate we excite dipole oscillations by 
abruptly shifting the center of the magnetic trap. With 
a single color lattice, the superfluid BEG oscillates freely, 
with the dipole frequency modified by the effective mass 
p7l.[l8j. Instead, adding an incommensurate bichromatic 
lattice, we expect the oscillations to be blocked by local- 
ization effects. Fig. 3(a) shows the center of mass motion 
after a shift in the magnetic trap of 6 /im, with a fixed 
number of atoms N = 1.5 x 10^, a fixed height of the 
primary lattice si — 10, and a variable height of the 
secondary lattice. At S2 = 0.1 (A/J = 2.6), the BEG 
oscillates with some damping. Increasing S2 the motion 
is strongly damped, until at S2 = 0.5 (A/J = 13), the 
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FIG. 3: (a) Measured dipole oscillations with increasing in- 
tensity of the secondary lattice. The height of the pri- 
mary lattice and atom number are fixed at si — 10 and 
— 1.5 X 10* respectively, (b) Measured dipole oscillations 
with fixed si = 10, A/J = 13 and two different atom number, 
Ai' = 1.5 X 10* and N = 2x 10^. The center of mass motion of 
the atoms was measured after 20 ms of expansion, fitting to 
the central peak. The dotted lines show damped sinusoidals 
fitted to the data with the frequency fixed. 



BEC does not move. Interestingly, the height of the sec- 
ondary lattice at which the condensate becomes localized 
is slightly greater than the non-interacting threshold for 
localization in the bichromatic lattice shown in Fig 1. 
The increase of this threshold could be due to the screen- 
ing effect of interactions shown in Fig. 2. 

We extended the investigations to a variable interac- 
tion energy effectuated by changing the number of atoms 
in the condensate. Fig. 3(b) shows the measured cen- 
ter of mass motion of the atoms for two different atom 
numbers TV = 1.5 x 10^ and = 2 x 10^, after an ini- 
tial harmonic trap shift of 6 /im, si — 10 and S2 = 0.5 
{A/J = 13). Transport is stopped only for the low num- 
ber of atoms, although damping is seen with the higher 
atom number. 

Localization due to disorder is not the only physical ef- 
fect which can block the motion of the atoms. In the case 
of an incommensurate bichromatic lattice the Bloch theo- 
rem cannot be applied, nevertheless the energy spectrum 
still shows the existence of energy bands. The band struc- 
ture for the single lattice is complicated by the emergence 
of 'mini-gaps' opening up almost everywhere across the 



spectrum [B"] . However the dominant modification to the 
first Brillouin zone of the primary lattice, for a weak ad- 
dition of the secondary lattice, are the extra energy gaps 
at kh and ki — kb, where ki = 27r/Ai gives the bound- 
ary of the first Brillouin zone of the primary lattice, and 
kh — 27r/Aft corresponds to the quasi-periodicity intro- 
duced on the larger length scale A;, = AiA2/(A2 — Ai) = 
4.38Ai from the beating between the two colors. This 
simplification to the energy spectrum is particularly true 
when interactions are introduced, since they can effec- 
tively screen the potential varying on length scales larger 
than the healing length, washing away the energy gaps at 
smaller k. As a consequence, in the presence of interac- 
tions one should carefully investigate the possible contri- 
bution of dynamical instability that in a sin gle lattice has 
been observed to block dipole oscillations [isl [T§ | when 
the quasi-momentum becomes greater than ~ O.bki. In 
the bichromatic lattice dynamical instability may occur 
at a small quasi-momentum ^ O.bki, corresponding to the 
beat periodicity, much smaller than the onset of instabil- 
ity with only the primary lattice. 

To better understand the contribution of the various 
physical effects instigating localization, the dipole oscilla- 
tions in our actual system were simulated using the time- 
dependent NPSE. Fig. 4(a) shows the center of mass 
motion using the same parameters as the measurements 
taken in Fig. 3(b); an initial magnetic trap shift of 6 
/im, si = 10, S2 — 0.5, and for both N = 1.5 x 10"* and 
N — 2.0 X 10^. The strongly damped center of mass mo- 
tion with low atom number and the increased movement 
with high atom number seen in the experimental mea- 
surements is confirmed by the solutions of the NPSE. 

From the momentum spectrum it is possible to distin- 
guish when dynamical instability is present. Fig. 4(b) 
shows the progression of the momentum spectra with 
N = 1.5 x 10^, and Fig. 4(c) shows the momentum 
spectra with = 2 x 10^. At t = 25 ms, the main 
components of the momentum spectra correspond to the 
peaks of the primary lattice at integer multiples of ±2?ifci, 
and peaks of the beat periodicity at integer multiples of 
±2hkh around the primary peaks. For both high and low 
atom number, at i = 50 ms extra momentum components 
begin to rapidly grow signifying the onset of dynamical 
instability [13]. We observe that the instability occurs 
when the quasimomentum becomes greater than 0.5fcb. 
By 75 ms there is a marked difference between the cases 
of high and low atom number. With N = 1.5x 10* atoms, 
the initial spectrum has been obscured by the additional 
momentum components. In contrast, with TV = 2 x 10^ 
atoms, the momentum spectrum retains much of the orig- 
inal spectrum structure. Intriguingly, the increased non- 
linearities inhibit the growth of the instability. 

Augmenting the interactions reduces the healing 
length of the condensate. In previous experiments, with 
a single lattice of spacing 0.4 /im and typically a healing 
length of ~ 0.3 ^m, the lattice spacing is not significantly 
greater than the healing length, and therefore measure- 
ments are largely indifferent to the atom number (isl . [l^ . 
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FIG. 4: (a) The center of mass motion from NPSE simulations 
of the dipole oscillations with si — 10, A/J — 13 and two 
different A*'. The dotted line shows the expected superfluid 
oscillation in a one color lattice with si = 10. (b) The corre- 
sponding momentum spectra in log scale at different times for 
= 1.5 X 10". The vertical arrows show the time at which 
each momentum spectrum is taken, (c) The momentum spec- 
tra in log scale at different times for N = 2 x 10^. 



However, in the case of the bichromatic lattice the large 
spacing of the beating (1.8/zm) allows interactions to ef- 
fectively smooth over the large scale beat periodicity. 
The growth of dynamical instability in our bichromatic 
lattice is governed by the competition between augmenta- 
tion by increased non-linearity and diminution by screen- 
ing of the beat periodicity [l^ . 

Consideration of the center of mass motion shown in 
Fig. 4(a) together with the momentum spectrum helps 
to unravel the different contributions to localization. The 
contaminated momenta at 75 ms for low atom number, 
shown in Fig. 4(b), is reflected in the complete blockage 
of the center of mass motion. However, most importantly, 
in the absence of dynamical instability, which is macro- 



scopically observed only after 50 ms, the movement of the 
atoms is strongly damped with respect to the superfluid 
case (see dotted line in Fig. 4(a)), suggesting that strong 
damping of the oscillations is not only due to dynamical 
instability originating from the beat periodicity of the 
bichromatic lattice. The presence of interactions and the 
possible screening of the minigaps at small k from the in- 
commensurate bichromatic lattice renders interpretation 
of the blocked motion at short times difficult. As already 
pointed out, with sufficiently strong interactions the en- 
ergy spectrum could be dominated by the periodicity at 
Af), and correspondingly the strongly damped motion at 
short times could be described by a high effective mass. 
Similarly to the screening behavior seen in Fig. 2, the 
contribution of the disorder will depend on the relative 
strength of interactions and the intensity of the disorder- 
ing potential. 

In conclusion, maintaining a minimal number of atoms 
we have observed a transition from oscillations to blocked 
motion with increasing intensity of the incommensurate 
bichromatic lattice. Simulations of our experimental pa- 
rameters using the NPSE show that the quasi-order in- 
herent in the quasi-periodic bichromatic lattice leads to 
the onset of dynamical instability, that contributes to 
the blocked motion only after a critical time. Screening 
of both localization due to the disorder and dynamical 
instability due to the beat periodicity was observed with 
strengthened interactions in the simulations. Increasing 
the number of atoms in the experiment, we observed a 
return of oscillating motion. 

This work shows that the exact choice of parameters 
is crucial to separate and isolate the effect of disorder- 
induced localization, or the non-trivial onset of dynami- 
cal instability in a bichromatic lattice. In future mea- 
surements one could focus on measurements at small 
quasimomentum to explore the transport behavior with- 
out the complication of dynamical instability, however in 
an interacting system the contribution of the disorder to 
blocked center of mass motion cannot be easily distin- 
guished from small amplitude oscillations resulting from 
an initial high effective mass when the energy spectrum is 
dominated by the beat periodicity. In the non-interacting 
limit, the bichromatic lattice is also a very promising tool 
to investigate Anderson-like localization that could be ac- 
cessed utilizing fermions or Feshbach resonances, which 
also provide the important possibility of tuning the in- 
teractions. 

This work has been funded by the EU Contracts 
No. HPRN-CT-2000-00125, IST-NoE-Phoremost, MIUR 
FIRB 2001, MIUR PRIN 2005 and Entc Cassa di 
Risparmio di Firenze. We thank G. Modugno for use- 
ful discussions and all the LENS Quantum Gases group. 



[1] B. Damski et al, Phys. Rev. Lett. 91, 080403 (2003); [2] R. Roth and K. Burnett, Phys. Rev. A 68, 023604 (2003). 



5 



[3] L. Sanchez-Palencia and L. Santos, Phys. Rev. A. 72 
053607 (2005). 

[4] V. W. Scarola and S. Das Sarma, Phys. Rev. A. 73 
041609(R). 

[5] R. B. Diener et at, Phys. Rev. A. 64 033416 (2001); 

[6] T. Schuhe et al, Phys. Rev. Lett. 95, 170411 (2005). T. 
Schulte et al, arXivlcond-mat/0609774 (2006). 

[7] J. E. Lye et al, Phys. Rev. Lett. 95, 070401 (2005). 

[8] D. Clement et al, Phys. Rev. Lett. 95, 170409 (2005). 

[9] C. Fort et al, Phys. Rev. Lett. 95, 170410 (2005). 
[10] M. Greiner et al, Nature 415, 39 (2002). 
[11] L. Fallani et o/., arXiv:cond-mat/0603655_ (2006) 



[12] P. Lugan et al, arXiv:cond-mat/0610389 (2006) 

[13] L. Fallani et al, Phys. Rev. Lett. 93 140406 (2004). 

[14] P. W. Anderson, Phys. Rev. 109, 1492 (1958). 

[15] R. T. Scalettar, G. G. Batrouni, and G. T. Zimanyi, 
Phys. Rev. Lett. 66, 3144 (1991). 

[16] L. Salasnich, Laser Phys. 12, 198 (2002); L. Salasnich, A. 
Parola, and L. Reatto', Phys. Rev. A 65, 043614 (2002). 

[17] This was confirmed in the experiment for si = 10 and a 
trap shift of 6 /^m. 

[18] F. S. Cataliotti et al, New J. Phys. 5 71 (2003) 

[19] A detailed analysis will be presented in a future theoret- 
ical work. 



